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Abstract 

The performance of frequency modulated continuous wave (FMCW) radar in tracking targets is presented and 
analysed. Obstacle detection, target tracking and radar target tracking performance models are developed and 
were used to investigate and to propose ways of improving the autonomous motion of unmanned surface vehicle 
(USV). Possible factors affecting the performance of FMCW radar in tracking targets are discussed and analysed. 



Introduction 

There have been growing interests in the use of fre- 
quency modulated continuous wave radar systems in 
USV applications such as autonomous navigation, tar- 
get ranging and direction (Chan & Judah 1998). Target 
range determination can be performed by means of 
well adapted frequency modulated continuous wave 
(FMCW) system. With FMCW, target range measure- 
ments can be made over large frequency bandwidth 
and FMCW system provides less accurate and unam- 
biguous range measurements (Reinhard & Schiek 
1997). In recent advancements, ultra-wideband fre- 
quency modulated continuous wave (UWBFMCW) 
radar systems have been developed to take advantage 
of microwave technology in providing simple solutions 
that can be adapted to suite the duration of modula- 
tion especially in tracking moving targets (Maaref 
et al. 2009). In recent developments, it has become a 
common practice to monitor coastlines with high fre- 
quency surface wave (HFSW) radar for sensing and 
monitoring ocean surface characteristics. The radar al- 
gorithms developed over the years were based on 
studies investigating the interaction between the ocean 
surface and radar signals. To achieve wider operating 
radar bandwidths for ocean scanning and target track- 
ing, linear frequency modulation waveforms (LFMW) 
are employed in such radar systems (Mahafza & 
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Elsherbeni 2004). Due to the high positioning accuracy 
and robustness of FMCW radar (Subramanian et al. 
2011), remote sensing of coastal lines using frequency 
surface wave radar was studied using FMCW as its 
primary wave form. This facilitated the extraction of 
target information on or close to the surface of the 
ocean using FMCW (Zhang et al. 2008). Obtaining 
target range autofocus characteristics required the 
critical consideration of the radar signal-to-noise ratio 
(SNR), radar dynamic range properties and geometric 
precision (Scheiblhofer et al. 2006). The linearity of 
radar autofocus characteristics can be achievable using 
phase-gradient algorithm and time-domain warping of 
de-chirped radar signal (Middleton et al. 2011). The 
non-coherent nature of radar systems has facilitated 
the development and integration low cost of radar 
sensors without the synchronisation of radar fre- 
quency or radar signal phase. Radar applications in au- 
tonomous navigation of unmanned surface vehicles 
(USVs) have introduced the classification of targets 
which are based on target range, cross-range, and tar- 
get velocity and radar power. This paper models and 
investigates the possible factors that affects the per- 
formance of FMCW radar in USV manoeuvring and 
target tracking. The description of radar waveform is 
provided. This is followed by obstacle manoeuvring 
and target tracking which provides insight to perform- 
ance monitoring of the radar in detecting and tracking 
obstacles. Finally the performance of the radar is eval- 
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uated in the presence of targets and noise while con- 
sidering the effects of radar sweep, power and the fre- 
quency of the waveform. 

FMCW radar wave form 

The choice of wave form and signal processing meth- 
odology in USV radar systems implementation was 
dependent on the type of mission assigned to the radar 
and the specific USV mission and function. The com- 
plexity and price associated to the development of USV 
systems constituted the major factor in deciding the 
type of radar waveform, hardware and software that 
were suitable for the specific functions that were 
assigned to the USV. Radar systems have the advantage 
of using continuous wave forms with or without fre- 
quency modulation. The frequency modulation of radar 
systems can be implemented either through analog 
techniques or digital methods. An important consider- 
ation in the modulation of radar frequencies was the 
radar range and Doppler frequencies as these had dir- 
ect associations to the choice of waveform frequency 
properties implemented in the radar system (Mahafza 
& Elsherbeni 2010). FMCW radar band pass signal x{t) 
can be presented as: 



x(t) = r(f) cos(2nf 0 t + O x (t)) 



(1) 



Where r{t) represents the radar signal amplitude 
modulation or envelop, O x (t) represents the radar signal 
phase modulation and f a represents the radar carrier 
frequency. The radar frequency modulation was mod- 
elled as: 



(2) 



The radar signal instantaneous frequency is modelled 



as: 



Jilt) =/0 +fm(t) 



(3) 



The radar signal x(t) may also be represented as an 
analytic signal forming the real part of the complex sig- 
nal if/{t) illustrated in equation (4) 

x{t) = Re{f(t)} = Re^rity®*^^- 1 } (4) 

The radar analytic signal is then defined as: 



(5) 



Where 

v (t) = rity^ 



(6) 



Implementing a Fourier transform to equation (5) 
while keeping the signal under the following conditions 
illustrated in equation (7) yields: 



ZX(a) 
0 



a > 0 
<y < 0 



(7) 



Where represents the Fourier transform of i//(t), 

gj = 2nf 0 and X(u) = is the Fourier transform of x(t). 
Implementing a step function to equation (7) yields: 



= 2U{a)X{a) 



(8) 



Where U(u) represents the step function of the radar 
signal in frequency domain. From the above models, it 
can be shown that if/(t) is represented as: 



f{t) = x{t) + jx(t) 



(9) 



Where x(t) represents the Hilbert transform of the 
radar signal x(t). The energy associated with radar sig- 
nal x(t) can be illustrated using Parseval's theorem 
(Mahafza & Elsherbeni 2010) as indicated below, 



1 



J x 2 {t)dt 



1 



(10) 



The exponential form of the radar signal can be mod- 
elled as (Zhang et al. 2008): 



i(t r ) = I o( j{«°tr+ane r ) 



<t r < 



T 

1 r 

Y 



(ii) 



Where I B represents the current supply to the radar, 
0) o represents the central radian frequency of the radar 
waveform and cj 0 = 2nf 0 , a represents the radar sweep 
frequency rate in Hz/s, t r represents the time variable 
and T r represents the radar sweep frequency interval. 
The radar frequency modulation sweep rate is the ratio 
of the radar frequency bandwidth B and sweep interval. 
Hence radar frequency rate is modelled as: 



B 



(12) 



The continuous variation in the transmitted frequency 
of FMCW radar at any given time with the difference in 
transmitted and received radar signal defines the radar 
beat frequency ft,. The radar beat frequency portrays the 
characteristic property and measure of target range R 
and it is modelled as (Chan & Judah 1998): 
cf b 
2/ m 



R 



(13) 
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Where c represents the speed of light in air, ft, repre- 
sents the radar beat frequency and f m = Af/T represents 
the rate of change of the transmitted frequency. The 
sweep frequency of the radar source is denoted by A/ 
and T denotes the time taken for each radar signal 
source sweep (Skolink 1981). The radar echo mixed with 
some portion of the transmitted signal received after r 
seconds produces the beat frequency of the radar. The 
radar echo is given as (Dorp & Groen 2010): 



Thus the instantaneous transmitted radar signal phase 
presented as (Dorp & Groen 2010): 



<p t (t) = ao + a\t + a.2 [ t 



To 
2 



(22) 



for0<t< T 0 

Having the beat signal phase (Dorp & Groen 2010): 



2r s 
c 



(14) 



Where r s represents the target slant range and c is the 
speed of light. The wave form transmitted by FMCW 
radar in compact form is given as: 



s(t) = Sj^ 



(15) 



<Pt{t)-4> t {t-T) 

To 



bo + bit + b 2 \ t 



for0<£< T 0 
With 



(23) 
(24) 



Where (p t represents the transmitted frequency phase 
with pulse width T, out power S 0 and bandwidth B hav- 
ing an instantaneous angular frequency of 



o) t (t) 



S4> t (t) 
St 



(16) 



The nonlinear transmitted frequency with instantan- 
eous up-chirp phase is modelled as: 



<P t (t)=2n{f 0 t + r -f 



(17) 



With the corresponding instantaneous frequency given as: 

(18) 



Similarly the radar instantaneous down-chirp phase 
and frequency is given as: 



fat- \* 



(19) 
(20) 



For -j<t<^-, where f 0 represents the radar centre 
frequency and fi = (2jtB)/T represents frequency modula- 
tion coefficient. The radar return echo with delay r 
modelled as s{t-r) is mixed with the transmitted radar 
signal to generate the radar beat waveform St,- 



S b (t)=s{t)s*(t-T) 



(21) 



bo = Cl\T + «2 T + 2fi2 r ~ 



b\ = 2fl2^ 



(25) 
(26) 



Obstacle manoeuvring and target tracking model 

The identification of obstacles and targets as USVs 
navigate on the ocean surface requires some form of 
target position deterministic algorithm which allows 
for target position and velocity estimation and deter- 
mination. To perform adequate autonomous motion 
operation, the USV is equipped and incorporated 
with obstacle and target tracking filters. The radar 
obstacle tracker provides the required optimum esti- 
mates of target position and velocity which are incor- 
porated in the dynamic model on which the target 
tracking filter is based on. The target tracking filter 
provides an accurate representation of the actual na- 
ture of obstacle motion path. The target tracking 
model is developed based on the assumption that the 
target manoeuvres at a constant velocity in which the 
straight line obstacle trajectory will be lost eventually 
should the obstacle manoeuvre in a different path. 
The target manoeuvring algorithm interacts with 
other USV autonomous motion deterministic algo- 
rithms in probabilistic models and it is based on Bar- 
Shalom's (Bar-Shalom & Cilang 1989; Bar-Shalom & 
Birmiwal 1982) interacting multiple model (IMM) al- 
gorithm. In the absence of motion characteristic 
properties of the obstacle, the obstacle is modelled as 
an object moving with constant velocity in a plane 
with the addition of process noise to incorporate 
slight changes in the target velocity. The resulting 
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dormant target model discretized over time T in the 
Cartesian coordinates is represented as: 



x(k + 1) = Fx{k) + G w{k) 



Where 



x 
x 

y 
y 

1 T 0 0 
0 10 0 
0 0 0 T 



(27) 



0 0 0 

w 2 

■ T/2 0 
1 0 
0 T/2 
. 0 0 



1 



(28) 

(29) 
(30) 

(31) 



With system characteristic process noise modelled as: 



£{w(/c)} = 0; E{w(k)w T (j)} = QS. 



(32) 



With system state estimate as X(0|0), system covariance 
as -P(0|0) and Qdkj detects whether target manoeuvre has 
occurred (Ramachandra 2000). 

FMCW radar performance evaluation setup 

To demonstrate the technique of target tracking with 
FMCW radar, the signal processing procedures integrated 
non-linear sweep frequency; bandwidth and power correc- 
tion techniques which were in Fourier transform for target 
range tracking. The FMCW upsweep frequency was linear- 
ized by the radar waveform generator to produce the phase 
modulated waveform modelled as: 



si(t) = a\Cos(p(t) 



And 



(p{t) = 2nf Q t + n — t 



(33) 



(34) 



Where <p(t) represents the radar phase at an interval t, 
fo represents the carrier frequency of the radar, B de- 
notes the radar sweep width and T m denotes the radar 
frequency pulse repetition period. The instantaneous 
transmitted waveform frequency is given as: 



F(t) 



1 d 

2n dt 



0M= fo 



B 



(35) 



With the radar frequency varying linearly at an inter- 
val of (0,T m ) between f 0 and fO + B. The FMCW radar 
waveform was split into two parts and passed into the 
radar mixer. The first part of the radar signal which is 
relatively small was used as a reference signal to facili- 
tate the detection of echo signals. The larger part of the 
radar waveform was passed to the circulator through to 
the radar antenna. The circulator provided multiport 
system that allowed electric signals to be propagated in 
clockwise direction as illustrated in Figure 1. The circu- 
lator facilitated the sharing of signals from the transmis- 
sion and receiving antennas. At the exit of the antenna, 
the radar waveform propagates out into the air towards 
an obstacle or target where it is reflected and returned 
by the receiving antenna. The received waveform then 
goes back into the circulator. An obstacle located at the 
distance r from the radar, generated an echo that is 
received by the radar waveform mixer in the form 
(Komarov & Smolskiy 2003) 



S2(t) = a2Cos(j){t-r) 



(36) 



Where r = 2r/c represents the radar echo propagation 
delay, c is the speed of light and «2 accounts for radar 




Signal Mixer 

Figure 1 Simplified FMCW radar circuitry. 
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waveform propagation losses, obstacle reflectivity and per- 
formance parameters of the FMCW. At the radar mixer, 
the echo signals and the radar reference signals are multi- 
plied. The process generates sum and difference waveform 
frequencies. The sum frequencies were usually in the 
order of twice the radar carrier frequency and the radar 
electronic circuitry cannot accommodate such high fre- 
quencies. Thus only the difference radar waveform fre- 
quencies were allowed to pass out of the radar mixer. The 
waveform that came out the radar mixer was modeled as: 



53(f) = «3COS[0(t)-0(f-T)] 

( r nB ' 

= a 3 cos\ 2nf beat t + 2nf 0 T- — t 

\ J m 



Where 



f, 



Bt 



beat 



' In 



dr 



(37) 



(38) 



And it represented the beat frequency of an obstacle 
echo. The radar pulse modulation frequency is denoted 
by f m = and dr = cl (2B) denotes the radar range 
resolution of the radar pulse (Brooker 2005). 

FMCW radar obstacle detection 

For autonomous motion of USV radar navigation appli- 
cations, obstacles are detected through the performance 
of Fourier analysis on the output waveform from the 
radar mixer. Implementing Fourier transform on the 
output signal with a rectangular window of integration 
over period of length T m generates the following power 
in the radar mixer output waveform: 



SechoififbeafB, T„ 



Pr 



sm[n(f-f beat )T m ] 

ntf ~f beat) 7 m 



(39) 



At the above model,/ represents the radar analysis fre- 
quency and P r = a\ represents the power of the echo 
from the obstacle. In using the standard radar equation 
to model the echo power, the echo power is modeled as: 



Pr=Pt 



oG t G r X 
(4tt)V 



(40) 



Where P t denotes the radar transmission power, a de- 
notes the radar cross section, A denotes the radar wave- 
length, r denotes the distance between the radar and the 
obstacle with G t and G r denoting the gains on the power 
scale of the transmitting and receiving antennas. In the 
Furuno radar used for the development of the USV, it has 
a single antenna and thus the gains (G t =G r =G) are equal. 
The use of non-uniform window of integration in the Fou- 
rier analysis of the output waveform from the radar mixer 
provided an effective measure in the signal analysis as it 



reduced side lobe effects in the spectral response of the 
radar signal. The implementation of a Hamming window 
of integration in the Fourier analysis provided means of 
detecting obstacles and targets that are within close range 
with large differences in echo level. The hamming window 
of integration for the Fourier analysis of the output signal 
from the mixer is modeled as: 



SechoifjfbeatjB, T m ) — P r 



(sm[n(f-f beat )T m ] 0.92 
V n(f-f beat )T m 1.08 



if ~f beat) T mSm[ n (f-fbeat) T » 



^-(f-fbeatfTl] 



(41) 



In implementing a rectangular window of integration 
at the output signal from the mixer, stronger echo at 13 dB 
will obscure weaker echo returned from an obstacle. Thus 
the detection performance of the FMCW radar is limited 
to the signal phase noise that was propagated through the 
circulator. The power spectrum of the signal phase noise 
(w/Hz) is modelled as: 



Sgen (f) = io 7 - 85 /- 305 



(42) 



Thus the FMCW beat signal from the mixer output with 
the inclusion of the signal phase noise is modeled as: 



s 3 (t) = a 3 cos[(p{t)-<p(t-T p ) + S<p(t)-S(p(t-T p 



(43) 



Where 8<p(t) represents the phase noise phase at time 
t, T p represents the difference in travel time on the para- 
sitic signal path at the circulator and the travel time on 
the reference path from the signal generator to the 
mixer. In the event that the parasitic signal path and ref- 
erence signal path are having the same travel time, their 
phase noise cancels out at that moment. The difference 
between the two process noises can be obtained using 
the transfer function model: 



H (f> T p) 



1-e 



(44) 



Where one of the noise processes is a time-delayed 
function of the other. The phase noise process limiting 
the radar performance is modeled as: 
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Sphase(f, Tp,L p ,) = P t 10 10 (\H(f, T p ) |) S generat e{f) 



P t lOW2[l-cos(2jifT p )]S generate (f) 



(45) 



The performance of the FMCW is limited by thermal 
noise having a power spectrum of 



J thermal 



(46) 



Where T denotes the temperature of the radar in degrees 
Kelvin, k B denotes Boltzmanns constant (1.38 x 10" JIC ), 
NF represents the noise factor in decibels. It is also an indica- 
tion of the noise factor decibel level increase of the thermal 
noise above theoretical lower limit. The signal-to noise ratio 
of the FMCW radar echo is modeled as (Bradley 2009a): 



SNR 



S echo {/ if beat? T m ) T m 
^p—\^ phase if i Tp^Lp^j -f- S thermal 



(47) 



The SNR is an indication of the power ratio of the 
sum of all limiting noise effects in the FMCW radar 
system having two parasitic phase noise paths. The de- 
tection of the radar echo is deemed reliable if the SNR is 
in the excess of 10 dB. 

FMCW radar performance prediction 

FMCW radar performance may be viewed as a two- step 
process. The first step takes into account the estimation of 
the statistical properties of the radar signal and noise fields. 
The characteristic property in this analysis is the signal to 
noise ratio (SNR). The second step converts the SNR to 
probability of detection using an appropriate mapping that 
is suitable to the radar type and mode of operation. If the 
SNR can be accurately determined, then the probability of 
detection is a number that lies on the interval {0, 1}. In 
practice, the SNR is not known with total certainty. 



Radar performance is predicted with aid of the stand- 
ard radar equation. Expressing the SNR ratio as a nor- 
mally distributed random variable x> the probability 
density function is modelled as (Bradley 2009b): 



h(x,fi, a) 



1 



exp 



2a 2 



(48) 



The SNR as a log-normal distribution y, its probability 
density function is modelled as: 



1 



10 1 



ln(10) y 



exp 



(10 log w {y)-(i) 

2cr 2 



(49) 



Where x and y have the relation, 



lOffl , X= 10 log w (Y) 



(50) 



The probability density functions represent the same 
information about uncertainty in SNR on two intensity 
scale. The scales are 0 < y < °o and - °° < x < °°. 

For optimisation of the radar having a signal with random 
phase and Rayleigh amplitude (Whalen 1971), the appropri- 
ate mapping SNR to probability of detection P D is given as 



P D {y,Pfa) = (Pfa) 1+ > 

With an inverse function of 



v(P D ,Pfa) 



H P fa) 

P D ln 2 {P D ) 



(51) 



(52) 



The inverse function has simple properties that pro- 
vide the platform of determining the probability density 



ShR = 15.4797 dB 



ShR = 15.4794 dB 



-100 



-150 




02 0.4 06 
freqjsiry (Nffe) 



02 04 0.6 0.8 
frequency (Mfc) 



Figure 2 Radar sweep analysis, (a) Target Echo @ Reduced Sweep Width, (b) Target echo @ Increased Sweep Width. 
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/be,, = 5.33333 MHz 



SNR = 23.£M4 dB 



! 





-150 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

0 12 3 4 5 6 , . /r 

frequency (Mfe) ^ ' 

(a) Sweep Width Investigation (b) Radar Beat Frequency Profile 

Figure 3 Radar transmission power, (a) Sweep Width Investigation, (b) Radar Beat Frequency Profile. 



5exp 



g(P D .P fa ,IA,0) 



10fa[-l+ln(p > )/ln(P 0 )] > 2 
r" 1 " ln(10) J l LO 



I H p f«) 



(rln(10)/« 2 (P D )[-l + \n(P fa )/]n{P D )]P D 



(53) 



function g{PD,Pfa,^,a) of the probability of detection P D 
in closed form. The solution is given as 

For P fa < P D < 1. a denotes the standard deviation of 
SNR (dB) and ^ denotes the mean of SNR (dB). 

As a result of the uncertainty in SNR, the probability of 
detection becomes a statistical distribution. The expected 
value and variability in the probability of detection are 
the useful characteristics of the statistical distribution. 



These characteristics are defined as: 
l 

(i Po = E[P D }= j P D g(P D ,P fa)f i,a) dP D (54) 

Pfa 

o 2 Pd =e[{P d -„ Pd ) 2 ] 

= J (Pd-Pd) 2 g(PD,Pfa,l*,o)dP D (55) 




0 1 2 3 4 5 6 Q0 01 Q2 03 04 05 06 0.7 

frequency (Mfc) frequency (Mfe) 

(a) Increased Sweep Frequency @Reduced Sweep Width (b) Reduced Sweep Frequency^ Increased Sweep Width 



Figure 4 Target echo analysis. 



Onunka and Nnadozie SpringerPlus 2013, 2:184 
http://www.springerplus.eom/content/2/1 /1 84 



Page 8 of 9 



u 



-fir 



□ 15.67 

□ 039 



SNR probability Atnaty 



04 

u 
u 

01 



- 



-5 0 5 10 15 20 25 30 0.0 
Figure 5 Probability of detection having lower signal standard deviation. 



p« probability dtruity 



■4pj ■ 0.J91 
<r. v .0.05« 



U 




10 



4 



Results 

The investigations on the FMCW radar revealed that the 
reduction of radar transmission power from 100 W to 10 
W reduced the obstacle echo power linearly slightly below 
the noise factor level (NF) in accordance to the radar 
transmission power as shown in Figure 2-a and also re- 
duced the frequency range of the obstacle echo frequency 
at the nominal sweep width of 5 MHz as illustrated in 
Figure 3. An increase in the sweep width of the radar 



resulted in an increase in the obstacle echo frequency. The 
radar beat frequency is attenuated at every second in a 
saw-tooth form shown in Figure 2-b. The further the obs- 
tacle, the less obstacle echo power was associated with the 
reduction in radar transmission power and radar sweep 
width as shown in Figure 4-a. Increasing the radar trans- 
mission power from 100 W to 1000 W and the radar 
sweep width increased the probability of detecting obsta- 
cles that are as far 5 km away as indicated in Figure 4-b. 
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Figure 6 Probability of detection having higher signal standard deviation. 
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In the analysis of the probability of detecting targets as 
shown in Figures 5 and 6, the distribution of the pro- 
bability of detection is peaked relatively small standard 
deviation when a is small. This implies that the predic- 
tion of probability of detection is useful and meaningful 
in this case. With a larger a, the distribution of the pro- 
bability of detection is much spread out and has much 
higher standard deviation. The prediction in this case is 
not useful as it is an indication of greater uncertainty in 
the signal. 

Conclusion 

The paper investigated the performance prediction 
model of an X-band 9410 MHz Furuno Radar emitting 
frequency modulation sweep at a 15 MHz short pulse 
and 5 MHz medium and long pulse repetition in obs- 
tacle detection and targeting tracking during maneuver- 
ing exercise of an USV. Obstacle range, radar cross 
section, the transmitting power of the radar and the gain 
of the radar antenna has noticeable effects on the per- 
formance of the radar in detecting and tracking obsta- 
cles. For effective performance of the radar to be 
achieved, it is proposed that the FMCW signal losses be 
reduced as much as possible and the time travel for the 
transmitted and received radar signal. The results 
showed that for a USV to be able to maneuver, detect 
and track targets effectively, the radar sweep width has 
to be increased to allow for narrower variability of radar 
signal and wider obstacle detection coverage and also 
increase in obstacle echo power and frequency. 
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